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Abstract 

On the one dimensional sphere, the support of the fundamental solution to the Schrodinger 
equation consists of finite points at the time t £ 2ttQ. The paper [Ka] generalized this fact 
to compact symmetric spaces. In this paper, we consider similar results on Zoll manifolds. 
We study the singularity for a solution to the equation using a functional calculus of the 
self-adjoint operator with integer eigenvalues. 



1 Introduction 

We consider the Schrodinger equation on Zoll manifolds. First we recall the setting of the problem. 
We study the following, 

(1.1) (idt + A M )u = 0. 

Here i — \j — 1 is the imaginary unit, t is the time variable and dt = Mr. M is an d-dimension com- 



df 

pact smooth Riemannian manifold without boundary. Am denotes the Laplace-Beltrami operator, 
i.e. in a local coordinate x — (x%, . . . , Xd) G M, 



d 

ij^i V det (5 4J ) V 



Here the matrix (g 13 (x)) > is the inverse of the metric (gij(x)) and det(gij) is the determinant 
of (gij(x)). Since M is compact, Am is a self-adjoint operator in L 2 (M), the space of the square 
integrable functions on M. By the spectral theory, we can write the solution operator to the 
equation (1.1) as e' itAjvf . We also assume M is the Zoll manifold, i.e. all geodesies has same 
minimal period. In a local coordinate, the geodesic flow is generated by the Hamiltonian vector 
field on T*M\{0}, 

(1.2) H = d i l-d x - dj ■ a 5 



were the function I is defined by l(x,£) — yX^f j=i 9 i: *( x )£i£j f° r x G £, ^ 0. Without loss of 
generality, we can assume the geodesic flow has minimal period 2-7T. Examples of Zoll manifolds are 
so called CROSSes (Compact Rank One Symmetric Spaces). For the background of Zoll manifolds, 
we refer [Be], jGulj . 

The following example is our starting point. 

Example 1. (^Ka , S. Doi^ Let M = S 1 ~ H/2nZ. We consider the following equation 
(1.3) 



{id t + d 2 x )G = 
G\ u=a = 6{x). 



Here S(x) denotes the Dirac measure at x = 0. Then the fundamental solution G(t,x) satisfies the 
following two statement 

(i) For fixed t/2n € Q, supp(G(i,x)) is finite points. 

(ii) For fixed t/2ir G R\Q, supp(G(£,x)) is the whole space S 1 
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Proof. First we prove the statement (i). By using the Fourier series, the solution can be written 
as follows 



(1.4) G{Lx) = — X/' 

77 feGZ 

Taking t = 27rn/m, n, to £ Z, and k = mj + I; j 6 Z, < I < m — 1. We have 

1 m— 1 

rn—l 

E —i27T — l 2 -\-ilx \ ^ imjx 
2ir^ 

i=o jez 
m— 1 m— 1 

■^-^ TO ' 

Z=0 j=0 

For final identity, we use the Poisson summation formula for 27r periodic functions. We define 

m— 1 m — 1 

(1.5) ff(n,m;i) = 1 £ e -^ 2 +^ = 1 £ 

TO ^— ' TO 

1=0 1=0 

We have the following representation of the solution 



m— l 

—i2iT — {mj+V) 2 +i{mj+V)x 

2n ■ 

j£Z 1=0 



m— 1 

(1.6) G(2Trn/m 1 x) = g(n,m;j)S(x — 2irj/m). 

3=0 

This imply (i). We prove (ii). We use the following symmetry 
Lemma 1.1. G(t,x) defined by (1-4) satisfies 

(1.7) G(t,x + 2t) = e i( ' x+ ^G(t,x) 

(1.8) G(t,x) =G(t,-x) 

as the sense of distributions. 

The proof of this lemma is a direct computation. Recall 

G = e ud *SeHi ( - sl ^\Hi(S 1 ) 

for any e > 0. Here H S (S 1 ), s € R are usual Sobolev spaces. Since S 1 is compact, there exists a 
point x such that near this point G(t,x) does not have regularity. So by (1.7), G(t,x) does 
not belong to H s near xo + 2t. Inductively we can prove similar singularity at xq + 2kt, k e Z. If 
t/2n — a e R\Q, these are dense in S 1 . Thus 

sing supp(G(27TQ!, x)) — S . 

□ 

Remark 1.2. In [Ka], Kakehi calculates g(n,m;j). He shows that for n,m € Z, to > 7 n,m: 
co-prime, 

g(n,m;j) ^ 0, Vj i/ to = 1, 3 mod 4, 

j g(n,m;j) = 0,forj:even ifm ^ 2modAj 

(1.9) [ 9{n 7 m;j) ^ 0, for j:odd, 

J g(n,m;j) ^0, for j:even 

< 1} m = V mod 4. 

[ g(n,m;j) = 0, for j: odd 

Moreover he generalized this example to compact symmetric spaces. 
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In this paper, we prove an analogous of this example. Before to say our result, we give a 
physical interpretation of this finiteness of the support. In quantum mechanics, the equation (1.3) 
represent the motion of the free particle. On S , the plain wave is given by 

This is the motion whose velocity is k. This means in S 1 the velocity is quantized. So at the time 
t, the particle can reach the point 

k ~ x = mod 2ir , k G Z. 
Taking t = 2irp/q, this can be written as follows 

P 

x = 2ir-k mod 2ir , k G Z. 

q 

These are finite points. By this naive interpretation, it is natural to expect same phenomenon 
for Zoll manifolds since \J— Am has asymptotically integer eigenvalues and whose geodesies are 
periodic. The main result of this article is the following. 

Theorem 1.3. Let M be a Zoll manifolds with period 2ir and xq G M . We consider the following 
equation 



(1.10) 



(id t + A M )u = 0. 
u\ t =o = <W 



Here S XQ is the Dirac measure on xq ■ The following two statement holds to this solution 

(i) For fixed t/2 € Q, sing supp(u(t, x)) is degenerate i.e. subset of finite d—1 dimensional 
sub-manifolds. 

(ii) For fixed t / 2it G R\Q, sing supp(u(£, x)) is the whole space M 

Remark 1.4. Moreover we compute the wave front set of the solution. We can also compute the 
wave front set to the equation with a first order perturbation. For more precise statement, see 
Corollary 3.4 and Theorem J^.l. 

We should remark some related results. In R n , the singularity propagates at infinite speed 
and the smoothing effect may occur. There are many results on the singularity of the solution to 
the Schrodingcr equation related to the smoothing effect. For example, if the classical trajectories 
are not trapped and the potential decays faster than the quadratic order, then the solution with 
compact support initial data is smooth at t > see e.g. |Naj . On the other hand to compact 
manifolds, there are few works to the study of the singularity except the work |Kaj . On R 1 , 
Yajima |Yaj proved fundamental solution is singular on the whole space time, to the super quadratic 
potential. Our result indicates similar result on Zoll manifolds. 



2 Preliminaries 

In this paper, we use the microlocal analysis on manifolds. We recall some relevant notations, see 
|Gr-Sj| , [Ho] . S m (T*M), rn e R is the space of functions a(x,£) on T*M, which are smooth in 
(x, £) and satisfy 

d^a(x,0 = 0((1 + |£ir^') ,(x,0 G T*M. 

We also write S-°°(T*M) = D m S m (T*M). If a G S m (T*M), we can define correspondent pseudo- 
differential operators, for short ^DO, by usual way. \E ,m (M) denotes the class of ^DO associated 
with this class. The smoothing operator is a operator which has smooth kernel. A operator in 
ty~°°(M) is the smoothing operator. It is useful to consider the operator as the operator on the 
half-density bundle. For the operator A G 1 $> m (M) and Riemannian density p 

Bu = piA(up-?), u G C°° (M, O') 
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defines B, the operator on the half-density bundle. We write this class of operator as \I> m (M; fl2,fl2). 
By this identity, we identify the operator on M with the operator on the half-density bundle. 

For A G y m (M;ni,ni), we have the symbol map a (A) G S m (T* M) / S m ~ 2 (T* M) and there 
exist a G S m (T*M), ax € S m - X (T*M) such that a (A) - a - a t G S m ~ 2 (T*M). We call 
a o = o~pri(A) the principal symbol of A and a± = a su b(A) the sub-principal symbol of A. Modifying 
above definition, ^t^[M\ Q,? , fis) denotes the class of ^>DO whose principal and sub-principal 
symbols are polyhomogeneous functions of order m and m — 1 respectively. \J — Am is a ^ DO 
G (M; Sis, fli) whose principal symbol is (Y^f J=1 {%)£,%£, j) 1 ^ 2 and sub-principal symbol is 
0. 

We also use the notation of the wave front set. If A G ^^ hg (M) has the principal symbol a, 
the characteristic set of A is given by 

CharA = {(x,£) G T*M\{0},a = 0} 

We define the wave front set of u G T>'(M) by 

WF(u) = n{CharA; A G *°^(M),Am G C°°(M)} 

where V(M) is the space of distributions on M. By definition, the wave front set is a conic subset 
of T*M\{0}. The wave front set is a refinement of the singular support, if H x (x, £) is the projection 
to x variable then 

n x (WF(u)) = sing supp(w). 
We use the following Propagation of singularity to the first order operator. 
Proposition 2.1. If A G 1 i ,1 (M;H,i has real-principal symbol. We have 
(2.1) WF(e" A «) = exp(tff a )WF(u),Vu G V'(M). 

Here exp(tH a )is the Hamilton flows for the symbol a = a~ pr i(A). 

We use a functional calculus to the self-adjoint operator with integer eigenvalues. By the 
following proposition, we can reduce y^Ajf to this type operator. This result is due to Colin de 
Verdiere, Duistermaat, Guillemin, and Weinstein-see |Hoj and |Gu2) for expositions of the spectral 
theory of the Laplace-Bertrami Operator on Zoll manifolds. 

Proposition 2.2. Let M be a Zoll manifolds with period 2tt. Then there exists a integer a and a 
self-adjoint operator E G ^~ hg (M;Cl^ ,Cl^) which commute with Am such that 



(2.2) Spec^ -A M + a/4 + E) C Z, 

For the proof, see Lemma 29.2.1 in [Hoj . a is so called the Maslov index, in our case this is the 
number of conjugate points in a closed geodesic counted with multiplicity. By this proposition, 
we should consider the self-adjoint operator with integer eigenvalues. The following special case of 
the functional calculus is essential to our results. 

Proposition 2.3. Let L be a self-adjoint operator on the Hilbert space H and f be a function 
Z — ¥ C. We assume Spec(L) C Z. Then we have the identity 

(2.3) f(L) = £ -L f(k)e^e^y L dy 



Proof. This is proved by 



2tt 
kez 



(2.4) J- r eW-^dy = 4,; 

27T Jo 

and the following identity 

(2.5) f(L) = £ f{k)P k ,e-^ L = e ~ lVkPk - 

kez kez 
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Here Sk,i = 1 if k = I and otherwise 5k,i = 0, Pk, k E Z are orthogonal projections to the space of 
eigenfunctions with eigenvalue k i.e. 



(2.6) 



Pk : projections to eigenfunctions with eigenvalue k, 

if k is an eigenvalue of L 
Pk : 0, if fc is not an eigenvalue of L. 



□ 



In the above identity, we shall change the sum and integral sign. We use self-adjoint operators 
(D x ) s on i 2 (S' 1 ) defined by the Fourier series, 



(2-7) (D x ) s f(x) = 5>) s / fc e' 



kx 
feez 



Here (fc) = (fc 2 + l) 1 / 2 , s e R and /(x) = £ feeZ /fce ifcx . This definition also can apply to the 
distribution on the 5 1 . We also define linear operators on H 

(2.8) (L) s = J2(k) s Pk. 

feez 

Since Pk are orthogonal projections, we get 

(D y ) N e- iyL = ^{l) N Pie~ ily 

(2-9) _ lEZ _ 

= ^(O^P ^ ^e"^ - (L) Jv e- <I ' i . 
(ez jez 

Proposition 2.4. Le£ L be a self-adjoint operator on the Hilbert space H and f be a function 
Z — > C. We assume Spec(L) C Z. T/ien we /iaue ifte identity 

(2.io) /(L) = (i) w ^ rj2( k y N f( k y kVe ~ tyLd y 

77 J ° feez 

ifJ2k£z(k)~ N f(k) etky * s absolute uniform convergent for some N G R. 
Proof. By Proposition 2.2 and (2.3), we have 

r-2-K 



/(£) = E f / I (r Ar w fcs w Jf E p ' e " i ' ! ' (i 
(2:m 

E^r [\k)- N f(ky k yJ2(i) N Pi^ Uy dy 

tew /7F •'0 



feez " u (ez 
Here Pfc are defined by (2.5). From (2.8), we get 

/(L) = (L)^ J] ^ 1^ {k)- N f{k)e-^ L dy. 
k% 2n J o 

Since J2kez( k )~ N f(k)e lky is absolute uniform convergent, we can change the sum and integral 

f{L) = r^ k r N tt k y kv z- ivL dy- 

2n Jo 



feez 

□ 
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Especially we have the identity 
(2.12) e - aL2 = (L) N ±- j 2 \D y )- N G{t,y)e-^ L dy 

for G(t,y) defined by (1.4) and N > 1. By (1.8), this can be written as follows 
e- UL ' = (L) N j-£(D y )- N G(t,y)(e^ L +e-^ L )dy 



(2.13) 

= IL) N 

2ir 

The following proposition is an abstract form of Theorem 1.3 (i). 



= • V t^1 \D y )- N G(t,y)co S (iyL)dy 



Proposition 2.5. Let L be a self-adjoint operator on the Hilbert space H. We assume Spec(L) C Z. 
Then for the time t — 2-jrn/m € 2ttQ, n,m <G Z, we have the following identity, 

m— 1 

(2.14) U(2Tm/m) = ^ g(n, m; j)V(2nj /m) 

3=0 

Here U(t) = e- ttL ' 2 , V(t) = e~ ttL and g(n,m,j) is defined by (1.5). 
Proof. By using the identity (2.12) and (1.6), we have 

r 2-K m—1 



U(2nn/m) = — / (D y )~ N V g(n,m;j)6(x - 2wj /m)(L) N e^dy 

= IT / £ 9{n,m;j)5{x - 2nj /m)(D v ) N e^ L dy 

= ^~ / 52 9(n,m;j)6(x - 2nj/m)e- iyL dy 

277 Jo j=i 



m— 1 

52 9(n,m;j)V(2nj/m). 

3=0 



Remark 2.6. Sy (%-13)> we can prove the similar result for V(t) — cos(itL). 



□ 



Remark 2.7. We can prove this proposition elementary from the following lemma which is proved 
by elementary linear algebra. 

Lemma 2.8. Let L and V(t) be as above and P; = J2jez Pmj+l-- There exists an invertible mxm 
complex matrix (a^), ay € C such that the following identity holds 

m— 1 

(2.15) p * = 52 OijV(2irj/m). 

3=0 

Later we use the wave front set of G(t,x) 
Proposition 2.9. Ift/2-n = a e i?\Q, we /iaue 

(2.16) WF(G(27ra, x)) = T*S\{Q}. 

for G(t,x) defined by (1.4). 

Proof. By example 1 (ii) , the singularity exists at x = 0. So we have (0,£o) € WF(G(27ra, y)) for 
some £o 7^ 0. By (1.8), We also get (0, — £o) € WF(G(27ra, y)). We apply similar argument used 
to prove Example 1 (ii). □ 
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3 Proof of Theorem 1.3 (i). 

In this section, we consider t/2-K G Q case. We prove the following theorem. 

Theorem 3.1. Let M be a Zoll manifold with period 2tt, P s G ^> phg (M;Cli be a self-adjoint 

operator and P = Am + P s - Then there exit self-adjoint operators Q\, Q2 G (M;Cli 

which is commute with Am and S G ^°(M; , such that at the time t = 2im/m G 2nQ, 
n,m G Z, the following identity holds 

m—l 

(3.1) U{2irn/m) = e~ iS W(2Trn/m) ^ g(n, m; j)V{2irj /m)e iS + Smoothing. 

3=0 

where U(t) = e ltp , V(t) = e~ ttQl and W{t) = e lt ® 2 . The principal symbols of these operator are 
given by o- prl {Q 1 )=a pri {^-A M ), o- pri (Q 2 ) = o- pri (P s + ay/-A M /2). Here 

1 f 2v 

(3.2) (T pri (P s ) = tt / <Tpri(P a )(exp(tHi))dt 

27T Jo 

where exp(tHi) is the Hamilton flow for the symbol I = o pr i{\J — Am). S satisfies 

e lS (A M + P s )e- lS = A M + P s + R 

for a self-adjoint operators P s G ^^ hg {M; fl^ , fis) which is commute with Am and R G \I/ _tx> (M; f2s , Q2 ). 
a denotes the Maslov index. The coefficient g(n,m;j) G C is defined by (1.5). 

Remark 3.2. e does not change the wave front set since it is invertible and it is a ^DO whose 
principal symbol is defined by e lapri ^ G S°(M). For V(t) and W(t), we can apply the propagation 
of singularity. 

Thus we have the propagation of the wave front set. 

Corollary 3.3. Let M and P = Am + P s as above. Then at the time t = 2nn/m, n,m G Z, 
m > 0, n,m: co-prime, we have 

(3.3) WF(U(2im/m)u) C exp(2im/mH q ) ^ exp(2nj/mHi)WF(u), Vu €V'(M). 

j;g(n,m;j)jtO 

Here cxp(tHi) and cxp(tH q ) are the Hamilton flows for the symbols I — (t(V-Am), and q = 
o{a.\J — Am/2 + P s ) respectively. 

The term j; g(n,m; j) ^ can be computed by (1.9). Similarly we get 

Corollaly 3.4. Let M and P = Am + P s as above. Then for the time t = 2irn/m, n,m G Z, 
m > 0, n,m: co-prime, the following identity holds 

(3.4) WF{U(2im/m)6 X0 )=exp(2Trn/mH q ) ^ exp(27rj/mfl-,)WF(J Xo ). 

Here 8 Xo is the Dirac measure at x G M and exp(tHi) and cxp(tH q ) are as above. 

This is a generalization of Theorem 1.3 (i). We prove Theorem 3.1 by reducing e ltp to the 
form to which we can apply Proposition 2.5. This is essentially well known but for readers we 
write the detail. First we will assume P s — 0. By Proposition 2.2, we take a self adjoint operator 
E G y~h g (M; n?,n?) such that 

(3.5) Spec( v /-A M + a/A + E)cZ 
We write 

(3.6) L = y/-A M + a/A + E. 
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Then we have 



L 2 = -A M + oV-Am/2 + E 

- phg V 



Here E £ ty® h (M;Cl' ,Cl') commutes with Am- Since E commute with Am, we have 



(3 7) e itA M _ e it(aV-A M /2+E) p -itL 2 

Next we shall consider the case P s ^ 0. First we assume that P s commutes with Am- Then by 
(3.7), we have 

(3 e itA u +Ps _ e it(uy/-A M /2+P s +E) e -itL 2 

Finally we consider the case that P s does not commute with Am- 

We claim that there exit self-adjoint 5 € $°(M;Sli ,Qi) and R € *"°°(M; fii, O^) such that 

(3.9) e lS (A M + P s )e^ s + R = A M + A- 

Here P s € ^^(M; fi^ ; fit?) is a self-adjoint operator which commutes with Am- The principal 
symbol of P s is given by (3.2). 

This is essentially due to [We] , see also [Ho] . Taking Q — P s — a \J— Am/2 — E, we will prove 
L 2 + Q is unitary equivalent to L 2 + B modulo smoothing operator where B G ^ , commutes with 
Am, is given below. We introduce the operator 

which is a self-adjoint ^DO. Taking an average, 

^ /-Su- 



it follows that 



Bi = — I Qtdt, 



o 

[L,B l ] = ±(e i2 * L Qe- i2 * L -Q). 
since 

at 

By (3.5), we have [L, B{\ = 0. Applying the Egolov's theorem the sub-principal symbol of L 2 + B\ 
is given by (3.2). We set B = B\. By similar calculus, we have 

(3.10) B-Q = 7T dt -rQsds = [IT, L] 

2tt Jo Jo ds 

where ^ 

T= — [ dt [ Q s ds. 
2^ Jo Jo 

Since L is elliptic, we can take its parametrix N. Then S\ = 1/2TN is a zero order self-adjoint 
^DO. By the calculus of ^DO, we have 

[5 1; i 2 ] = \{[T,L 2 ]N + T[N 1 L 2 }) = [T,L] + R 

where Rq £ 'J' -1 is a self-adjoint operator. We set S = Si and consider ie zS (L 2 + Q)e~ lS . 

If A € \t m , then ie Ae~ lS is a ^DO whose symbol is the asymptotic sum of the following 
formal scries 

oo 

]T((ad iSyA/j\ 



where (ad iS)A = [iS, A}. And by (3.9), We have 

ie lS {L 2 + Q)e- iS = L 2 + B + Ri 
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where R\ £ "J is a self-adjoint operator. Applying similar argument to R±, we have B 2 € 'J and 
S2 € satisfying 

ie iS (L 2 + Q)e~ l5 = P 2 + B + R 2 

where B = B\ + B2, S = Si + S2 and i?2 € W -1 are self-adjoint operators. Continuing this 
argument and taking asymptotic sums, we have desired B and S. This means (3.9). Thus we have 
the identity 

(3.11) id t + A M + P S = e lS (id t + A M + P s ) e - lS + R. 

We shall give the solution operator to this equation. We write the solution operator as follows 

and inserting this to the right hand side of (3.11), we have the following equation to the operator 
A(t), 



(3.12) 



((id t + A AI +P s + e- lS Re lS )A(t) = -e^ 5 Re lS e «(A M +P.) 
|A(0) = 0. 

This can be solved by Duhamel's principle, we have 



(3.13) Mt)=i e^- s ^' +p ^ lSRe '^e- lS Re lS e ls ^' +p ^ds. 



Since e lS Re zS is a smoothing operator, A(t) is a smoothing operator. We have the identity 

(3.14) 6 «(Am+p s ) = e -iS e it(A M +P s ) e iS + Aif y 

Since P s is commute with Am, by the identity and (3.8), we have proved the following 

Proposition 3.5. Let M be a Zoll manifold with period 2tt and P = Am + P s as above. Then 
there exit self-adjoint operators S, E £ ^ (M; Qz , fi? ) such that 

(3.15) e it(A M +p 3 ) = ^s^^K^n+Ps+E)^^ + Smoothing 
Here L is defined by (3.6) 



We apply Proposition 2.5 to (3.15). Taking Q\ = L and Q 2 = ct\J— Am/2 + P s + E, we have 
Theorem 3.1. 

Remark 3.6. If the manifolds are odd dimension rank 1 symmetric spaces, we can easily prove 
support theorem by our identity although Kakehi proved support theorem for more general compact 
symmetric spaces. 

We outline the argument. Essentially we can assume the manifolds are odd dimensional spheres. 
The spectrum of the Laplace-Beltrami operator is the following, see e.g. |Taj . 

(3.16) specA Sti = {-k(k+(d- l)/2);fc e N}. 

We have 



A s , = A sd - (d - l) 2 /4 + (d - l) 2 /4 = -^-A sd + (d- l) 2 /4 + (d- l) 2 /4 
By (3.16), we know 

spec V-Asd + [d - l) 2 /4 c Z 

So we can apply the Remark 2.5 to this operator. Using the Huygens' principle, see |La-Phj . |Taj . 
we have the following 

(3.17) suppte 2 ™"/™^^,) = J2 d(27rj/m)supp(6 X0 ),WueV'(M). 

Here S Xo is the Dirac measure at xo £ S d , d(t) is defined by 

d(t){x} = {ye S d ;dist S "(y,x) =t}, t G R,.t e S n , 
where dist^d is the distance of the Sphere. 
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4 Proof of Theorem 1.3 (ii). 

We shall consider t/2-K £ R\Q case. We prove the following a generalization of Theorem 1.3 (ii). 

Theorem 4.1. Let M be a Zoll manifolds with period 2tt, P s g 4^ (M; f^ 5 , fi^) be a self-adjoint 
operator and P = Am + P s - If t /2ft = a G R\Q, t/ie following identity holds 

(4.1) WF([/(2to)4 ) = exp(27rai/ 9 ) U < T < 2 , cx P (tH p )WF(5 Xo ) 

Here 5 XQ is the Dirac measure at xq £ M, U(t)S X0 — e ltp 8 XQ , exp(tHi) and exp(tH q ) are the 
Hamilton flows for the symbols I = o~ pr i(y/— Ajv/ ), and q — cr pr i(P s ) respectively. cr pr i(P s ) is 
defined by (3.2). 

Remark 4.2. We can also prove similar identity to WF^/ 2 - We notice the following 

(4.2) WF(U(t)u) C exp{tH v ) U < T < 2w exp{ T H p )WF{u) 
by (2.12) and Proposition 3.5. 

Remark 4.3. Since M is complete, if o~ pr i(P s ) = 0, then we have 

sing supp(t/(27ra)(S a;Q ) = M 

from Theorem 4-1- This implies Theorem 1.3 (ii). 

By Proposition 3.5, we study the singularity of e~ ltL S Xo here L is defined by (3.6). We fix 
the coordinates (x, £) of the cotangent bundle where x is the normal coordinate around xq = 
and £ is the corespondent coordinate. We study e~ ttL S Xo near xo. We recall that e~ ltL can be 
approximated, modulo smooth functions, by the Fourier integral operator, see |Gr-Sj| , [Ho| 

(4.3) e ltL u=(2Tr)- d / 2 [ e i ^ t ' x '^- y ^a{t,x,Ou(y)dyd^ mod C°°{M). 



From now on, we omit the term mod C°°(M). Here the phase $(<, x, £) and the amplitude a(t, x, £) 
are chosen to satisfy 



(4.4) 



(id t + L)e^- x ^a(t,x,0 = 0((1 + 
e»(°.*.0a(O,a;,£) = (2ft)- d / 2 e ix <. 



These functions can be taken if <!>(£, x,£) satisfies the following Eikonal equation 
(4.5) 



d t $ + l(x,d x $) = 
$(0,a:,£) = af-£ 



here the symbol l(x,£) is the principal symbol to the operator L and we also inductively choose 

the amplitude a(t, x, £) = do + ffli H by solving transport equations. For small i, these equations 

have solutions and we get $(t, x, £), a real C°° function homogeneous of degree one in £ and 

a(t, x, £) = a + aH e Sphg(M) where € S ,_J (M) is homogeneous of degree — j in £. In (4.3) 

we take u = S Xo , and get 

e ltL S Xa = J e i *<t'*'*'>a{t,x,S)dZ =: A(t,x). 

Here the integral is defined by the oscillatory integral. Taking a point x\ =t xq near xq, we study 
A(t,x) in some neighborhood of x%. Change coordinates to polar coordinates (r, 6), (A,w) to x, £ 
respectively. We have 

A(t,rd)= I e I * ( *' re ' Aw) a(t,r0,Aw)A d - 1 dAdw 

JQ JSU-l) 
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We write x\ — r\6\. Here n > is sufficient small. We shall apply the stationary phase method 
to the lj integral. Since $ satisfies (4.5), we have 

On normal coordinate, the metric satisfies 

(4.6) g i '(x)=5 ij + 0(\x\). 
We take \x\ = r the same order as t. Then we get 

(4.7) x,0=x-€- t\£\ + O(t 2 = X(r9 ■ lj - t + 0(t 2 )) 

Since 9 ■ lj take its maximum or minimum at lj = ±9 respectively, lj = ±9 is critical values of 
9 ■ uj. Moreover they are only critical points and non-degenerate i.e. the Hessian is not degenerate. 
So we can assume the critical points of &(t,r9, Xlj) is non-degenerate. By the stationary phase 
method, we obtain 

A(t, r9) = J2 e iA *^ r %(i, rO, \)\W 2 dX 
The phase function is defined by 

$j(t,r9) = $(t,r6,LJj) 
where ljj are finite critical points to <&. The amplitude is 

ajit, r6, A) = fijo + Oji + ■ ■ ■ G S° hg (U) 

considering t and 9 as C°° parameter and U is some neighborhood of n . We have $j = (r — t) + 
0(rt) or $ = — (r + t) + 0(rt) from (4.7). We also assume d r $j ^ by changing r ~ t smaller. If 
there exists a point r 2 6 l 2 such that &j(t, Ti9i) ^ then by using the identity 

-Q^dxe"*' = e ixi * 
i J 

and integration by parts, we get the integral is C°° near this point. On the other hand, by the 
propagation of singularity, A(t, x) has only singularity at r = t, —t. So essential parts of the integral 
are the phase $ 3 s which are the following form 

$+ = (r - t)$+@or@$~ =-(r + t)$j. 

Since d r $j ^ 0, we have $^ ^ 0. Near n, we get 

A(t,r9)=Y] / e lX ^ t ^t^d 3 (t,r9,X)X^/ 2 dX 

pOO 



By changing the variable as A' = A3>^(£, r9) 

poo 

A(t,r9)= / e iX '( r -Va+{t,r6,X')X' 
Jo 

e- iX ^a-(t,re,X')X' {d - 1)/2 dX' 
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where a± = ^ . a,j(t,r6, \<&j)&j . By the propagation of singularity, we know a± is elliptic. If 
t > 0, since (r + t) ^ the second term is C°° , we can write this integrals as follows 

/>oo 

4(t,r0)= / e iA '( r -*)a + (t ) re,A / )A /(< '" 1)/2 dA / 
io 

5* x(r -^a' + {z, r6, \')\ /{d - 1)/2 6(z - t)dzd\' 



A(t,r9) = I e- iX ^a-(t,re,X')X' id - 1)/2 d\' 



oo poo 



'0 Joo 

If t < 0, the first term is C°° and we have 



/0 /-oo 
/ e 4A '( r - z )a'_(z,r0,A')|A'| ( ^ 1)/2 (5(z- |t|)tfadA' 
-oo J—oo 



Here a'_ = a_(z,r#, |A'|). Thus we have 

Proposition 4.4. Let M 6e a smooth Riemannian manifold, xq G M and L is a self-adjoint ^ DO 
satisfies a pr i(L) = o pr i{\J — Am). W^e tafce a;i 7^ mo sufficiently near xq and the normal coordinates 
r6 around x = 0. T/ien £/iere exists U, a neighborhood of r\, and a'±(z,r6, A') € S^ hg (U), elliptic 
symbol with C°° parameter 9 such that for small t > 0, we have the following identity near X\ 
modulo smooth functions 



(4.8) e ±ltL 5 Xa {r6) = I 



pOO pOO 

/ / e iy ^d' + (z,r9,\')\' {d - 1)/2 6(z-t)dzd\' 

JO J -00 

/0 pOO 
/ e lX ' {r - z) a'_{z, r9, X')\X'\ {d - 1)/2 S(z - t)dzdX' 
-00 J —00 

Remark 4.5. We can get this fact more directory by solving (4-5). For short, we use the propa- 
gation of singularity. 

Now we study the singularity of e~ ztL ' ' 5 Xo near x\. We take a smooth cut-off function p{r9) 
near (n, #1), i.e. p = 1 on a small neighborhood of (n, #i) and p is identically zero out side another 
neighborhood of (ri, Then by the identity (2.13) and the propagations of singularity, we have 
the following identity for N > 1 

(4.9) P e-* tL2 S X0 ee p(L) N p p(y){(D y )- N G(t, y)}{(e^ p + e^y p )}d y S Xo 

Jo 

where p(y) G C°°([0,7r])) is a cut-off function near r . By Proposition 4.4, the right hand side can 
be written as follows in the normal polar coordinate near xo- 

(e iyP + e-^ p )5 X0 ee — / / e iX '^b(z, r6, X')5(z - y )dzdX'} 

^ 7r J — oo J — oo 

Here we take 

in x ui ... \2TTd' + {z 1 re,X l )X l[d - 1)/2 A'>0 

(4 ' 10) 6( ^ M) ^2^Lm')|a'|^ 2 A'<0. 

Since a!± is elliptic with respect to A' , b is elliptic. The right hand side of (4.9) becomes 

pTT i /'OO /*oo 

p(L) N p(y){(D y )- N G(t,y)}— / e^"*)^, r0, A')*(* - i/)<fedA'd» 

J ^ J-co J-co 

1 /*00 /*00 /*7T 

= p(L) N — / e* A >- z >6(z,rM') / p(y){(D y )- N G(t,y)}S(z - y)dydzd\' 

Z7r J-oo J-oo JO 

i /*00 pOO 

= P(L) N - / e iA '('-^K^^,A'){p(y)( J D !/ )- JV G(t )2/ )}| ?/= ,^dA' 

^ 7r J — oo J—oo 

Here the last integral is *DO with C°° parameter 9. We write this operator as B{9) and get 
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Proposition 4.6. Let M be a Zoll manifold with period 2n, x <G M and L is a self-adjoint ^>DO 
defined by (3.5). We take X\ 7^ x sufficiently near x and the normal polar coordinates r9 around 
xq = 0. Then there exists U, a neighborhood of r\, and b(z,r$, A') € S^ 1 ^ 2 ^), elliptic symbol 
with C°° parameter 9 such that for small t > and n > 1, we have the following identity near x\ 
modulo smooth functions 

(4.11) pe- itL 'S X0 = (L) N B(9)(p(y)(D y )- N G(t, y)) 

This proposition reduce the singularity of e~ ltL 5 Xo to G(t, y). If t/2n = a € R\Q, by Proposi- 
tion 2.9, G(y) = p(y)(D y )~ N G(2na,y) is singular at r\ to both directions. Since B(t,9) is elliptic 
'J'DO with C°° parameter 9, this preserve the singularity to r-variable. We have 

(r 1 9 1 ,±9 1 )eWF(e- ltL2 S X0 ). 

We also know, in normal polar coordinate, 

(4.12) (ri0i,±0i) = exp(±riir,)(O,±0i). 
Since we can take x\ arbitrary near x , we have 

(4.13) (O.iOiJeWFfe-^^). 

Next we prove the existence of the similar singularity at every points using (1.7). As in the previous 
argument, we take smooth cut-off functions p(r6) near n x (exp((n + An a) Hi) (0,9 1)) and p(y) near 
ri + Ana We have 

pe- i2 ™ L2 6 X0 = p(L) N f p(y){(D y )- N G(2na,y)}(e^ L + e-^ L )dy5 Xa 
Jo 

= p(L) N ^ p(y){(D y )- N e^ a +y^G(2na, y + Ana)}(e lyL + e~ iyL )b Xa dy. 
Jo 

Change the variable y = y' — Ana. 

= p(L) N f p(y' - Ana){(D y )- N e l(y '- 27T ^G(2na,y')}(e l(y '- 4 ^ L + e-^ y '- 47Ta)L )S Xo dy' 
Jo 

= pe l4naL (L) N [\p(y' - Ana)(D y )- N e^'- 2 ^G(2na, y')}e iy ' L S Xa dy' 
Jo 

+ pe-^ aL (L) N l\p(y' - Ana)(D y )- N e^'- 2 ^G(2na, y 1 )^' L 5 Xa dy' . 
Jo 

We write this as follows 

(4.14) pe- l2lTaL2 5 X0 = pe l4mL l + pe- ;AlxaL II. 

Now G(y') = p(y' - Ana)(D y )^ N e l(y ' ^ 27Ta '> G(2na, y') is singular to both directions at y' = n 
and applying previous argument. We write (I) as B\(9)G(y'). Since B\(9) a ^DO elliptic to the 
positive direction r > 0, we have 

(ri0i,0i) e WF(I) 
Taking n sufficient small, we also assume 

— n — 8na ^ n + 47ra, mod 27T 
By shrinking the support of p, we have 

exp(-87rai?i)(ri0i, -ft) £ WF(I). 
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By the propagation of the singularity, we get 

, . expf^affOtnfbfiJeWF^ 4 ™^) 

( 4 - 15 ) -a r 

exp(-4 7 ra J ff ( )(r-i6'i, -ft) £ WF(pe li * aL 1). 

Similarly we obtain 

exp(Wff ; )Mi,#i) t WF(pe- i47rai II) 
exp(-47ro!^)(ri0i, -6>i) e WF(p e - i47raL II). 

(4.12), (4.13), (4.14) and (4.15) mean 

exp((n ±47ra)Hj)(0,±«i) G WF( e - 42 ™ L X)- 
Since n > is any small constant we have 

exp((±47ra)ff,)(O,±0i) e WF(e- i27rai X)- 

Inductively we can prove similar singularity to the points exp((47rafc)iJ;)(0, ±#1), Vfc e Z, V0i. 
Since Airak, € Z are dense in S 1 , we get 

(4.17) U < T <2. exp(Tflp)WF(<5 X0 ) C WF( e - l2 ™ i2 6 Xo ) 

Here we use periodicity, exp(-tHi) — cxp((27r — t)H{). By Proposition 3.5 and (4.2), this implies 
(4.1). 
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